SIViP (2015) 9:1975-1982 
DOI 10.1007/s11760-014-0692-x 


ORIGINAL PAPER 


C) CrossMark 


A fast line-scanning-based detection algorithm for real-time SAR 


ship detection 


Xiaolong Wang - Cuixia Chen 


Received: 20 July 2013 / Revised: 27 August 2014 / Accepted: 27 August 2014 / Published online: 11 September 2014 


© Springer-Verlag London 2014 


Abstract Synthetic aperture radar (SAR) provides a pow- 
erful surveillance capability allowing the observation of tar- 
get, independently from weather effects and from the day 
and night cycle. Unfortunately, the automatic interpretation 
of SAR images is often complicated and time consuming. In 
support of real-time vessel monitoring, a fast line-scanning 
detector designed for detecting ships from SAR imagery is 
proposed in this paper. The detector does not require any 
prior knowledge about ships and background observation. 
It uses a novel local gray-level gathering degree algorithm 
to detect potential targets and then a complementary filtering 
scheme to reject false alarms. The performance analysis over 
real SAR images confirms that the proposed detector works 
well in various circumstances with high detection rate, fast 
detection speed and perfect shape preservation. 


Keywords Local gray-level gathering degree (LGGD) - 
Real-time - Ship detection - Synthetic aperture radar (SAR) 


1 Introduction 


Ship detection is of great significance for various mandates 
associated with marine monitoring, such as traffic surveil- 
lance, military reconnaissance, fishery monitoring, as well 
as vessel search and rescue [1]. With the increasing volume 
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of image data that are collected from air- and space-borne 
SAR systems, a large amount of investigations have proven 
that SAR can be used for ship detection purpose. During the 
last two decades, much effort has been devoted in develop- 
ing algorithms, establishing infrastructures and even build- 
ing automatic target recognition (ATR) systems for process- 
ing marine SAR images [2-8]. The distributed constant false 
alarm rate (CFAR) framework is the most widely accepted 
conceptual model. CFAR-like detectors, as the widely used 
ship detection approach on which several current operational 
systems rely [5], involve the parameter estimations of ships 
and local clutters. Also, the threshold setting is essential so 
that a constant false alarm probability is guaranteed for all 
values of unknown clutter parameters. The strong depen- 
dence of the CFAR schemes on prior knowledge about ships 
and background observation limits their application in auto- 
matic detection. Besides, due to differences in environmental 
conditions such as changes in clutter edge, multiple targets 
or jamming, the target detection is often corrupted. 

The development of maritime ATR system aimed at imple- 
menting a simple, yet effective means of rapidly alerting and 
cueing the supervisors to regions within high-resolution SAR 
imagery that could contain targets of interest [9,10]. Auto- 
matic detection is the key step of the maritime SAR ATR 
system. Although great effort has been paid, fast and robust 
ship detection still remains a challenge. Currently, the devel- 
oped algorithms mostly emphasize the detection rate but not 
the time, without regard to construct a practical system, just 
as some CFAR-like methods. Actually, real-time processing 
is usually a requirement for the operational ship detection 
systems. In recent years, Duman et. al. used region covari- 
ance (RC) and codifference algorithms to solve SAR ATR 
target detection problem [11,12]. The proposed method is 
shown to deliver high detection accuracies and low false 
alarm rates. However, its computational complexity is still 
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high due to RC algorithm, codifference algorithm and train- 
ing and testing steps with support vector machines (SVMs), 
although directional filters have been introduced to decrease 
the search space. Computational complexity is a major prob- 
lem in developing a fast detection algorithm. High compu- 
tational complexity always causes a large computation cost, 
hence, affecting the system’s implemental efficiency. There- 
fore, it has become an urgent problem to improve the detec- 
tion speed on the condition of high detection rate. None but 
simple and fast ship detection algorithm can afford to meet 
the expanding requirements of intelligence, surveillance, and 
reconnaissance. 

Motivated by the characteristics of spatial intensity distri- 
bution as observed in any SAR images, a “local gray-level 
gathering degree (LGGD)” algorithm for ship detection is 
developed. Based on this, a fast line-scanning detector is pro- 
posed in support of real-time vessel monitoring. The detector 
constituted by a linear LGGD algorithm and a postpositive 
complementary filtering scheme, is named S-LGGD detec- 
tor. The merit of a postpositive filtering scheme rather than a 
prepositive one (which is generally adopted by most detec- 
tors) is the latter always reduces target information during 
noise suppressing. Performances of the detector have been 
investigated on several real SAR images, and moreover, a 
practical real-time ship detection system for an X-band air- 
borne SAR system has also been established. Experiments 
show that the proposed detector works well in various cir- 
cumstances with high detection rate, fast detection speed and 
perfect shape preservation. 


2 Algorithms and methods 


Cell detection is the first canonical problem to be treated in 
developing a framework for automatic exploitation of ship 
information from SAR imagery. Many systems rely on the 
CFAR approach which interprets the difference between the 
target and its surrounding ocean clutter as a difference of 
intensity levels. However, the validity of the CFAR approach 
strictly depends on assumptions, hypotheses and even prior 
knowledge about ships and background observation, which 
increase its computational complexity and put its applica- 
tions into a dilemma of the choice between performance and 
speed. To solve the problem, the new detection algorithm 
must exhibit low computational complexity and a reasonable 
false alarm rate. 


2.1 Description of LGGD algorithm 
The intensity of each resolution cell in SAR imagery denotes 
the radar scattering characteristics of arandom ground object. 


Cells with different intensities and spatial distributions often 
represent the various ground objects. Suppose that a cell 
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Ci (i, j) with intensity J has the same observed intensity with 
another cell C2(i+ Ai, j + Aj), they are coherent and there- 
fore the ground objects represented by the two cells should be 
the correlative objects. Accordingly, a gathering relationship 
of the cell pair (C1, C2) at position (i, j) can be expressed 
by 


G = 1/ (exp(Ai) + exp(Aj)), (Ai)? + (Aj)? #0), 0) 


where Ai and Aj represent the spatial intervals of the cells in 
two directions (range and azimuth), respectively. For all cell 
pairs in relation to the cell Ci (i, j), a gray-level gathering 
relationship of the intensity 7 at (i, j), namely gray-level 
gathering degree, can be computed as 


M N 
GrG = $, $, 1/(@xp(Ai) + exp(Aj)), 
Ai=1 Aj=1 
(CG + Ai, J+ A) =D, (2) 


where M and N represent the spatial dimensions of the SAR 
image known as azimuth and range, respectively. 

The G; describes the spatial distribution of cells and inten- 
sities, and hence, it is able to highlight the changes of the 
spatial information. Unfortunately, its computational cost is 
increased with the image dimensions. To decrease compu- 
tational cost, a LGGD is adopted by limiting a calculation 
range (window). It is obvious that the LGGD algorithm is 
related not only with the number of matched cells, but also 
with their spatial distributions. If the cells with gray-level 7 
come from the same surface of a ship target, they are gather- 
ing and hence Gy is large, while if they come from the ocean 
clutter noise, they are discrete and then Gzy is small. Thus, 
ships should be detected if an appropriate threshold Gih is 
chosen. 


2.2 Design of S-LGGD detector 


The LGGD algorithm does not require any prior knowl- 
edge about ships and background observation. The speed 
and intensity sensitivity are its outstanding merits. With these 
advantages, a novel ship detector is proposed in considera- 
tion of the line-by-line data downlink and processing mode 
for most real-time SAR systems (especially for airborne SAR 
systems). The block diagram of the developed detector is 
shown in Fig. 1. A linear sliding window is designed for per- 
forming a line-scanning detection. The window constituted 
by guard windows and reference windows includes 2N + 2L 
range samples (range cells) corresponding to 2N reference 
cells and 2L guard cells surrounding the test cell. Due to 
lack of ability for distinguishing bright cells from dark cells, 
the LGGD algorithm cannot distinguish from false alarms 
corresponding to speckle noises, sea clutters and even dark 
areas (such as low-wind areas, oil spills or rain cells). As a 
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Fig. 1 Structure of the proposed detector (S-LGGD detector). 


countermeasure, a complementary filtering scheme of gray- 
level discrimination and isolated point removal is designed to 
reject the false alarms. The detection process can be divided 
into three steps: 


(1) LGGD detection: For each inputted line, a sliding detec- 
tion is executed by computing the LGGD value of each 
cell with its 2N reference cells. Owing to the higher 
intensities of ship cells, there are often no matched cells 
in the sliding window for them, that is, a zero value of 
G; is universal for ship cells. Therefore, potential target 
cells corresponding to the zero value of G z are presented. 

(2) Gray-level discrimination: For each potential target cell, 
an auxiliary gray-level discrimination based on its image 
intensity is made, and then, target cells are differentiated 
from false alarms. 

(3) Targets confirmation: For the remaining target cells, 
through a supplementary noise elimination, the isolated 
cells are removed and the ship cells are confirmed. 


The window size is vital to the real-time LGGD detention. 
Too large size is disadvantageous to calculation efficiency, 
while too small size is not enough to calculation stability. In 
order to obtain a reasonable size for automatic detection, we 
have investigated relationships of the LGGD value with the 
increasing window size by statistical analysis method. For 
the guard window, comparisons derived from several SAR 
images under different background noises have shown that: 
Five pixels window size (L = 5) is enough to protect the test 
cell, and more pixels could hardly contribute to a different 
result. For the reference window, a series of theoretical max- 
imum LGGD values with the increased cells are calculated 
on the assumption that their internal cells are all coherent. 
Due to the exponential operator in the LGGD algorithm, the 
calculated values tend to be invariable when reference cells 
are over 20 pixels. Thus, the window sizes of 30—40 pixels 
(2N =20-30 and 2L = 10) are suggested in our scheme, 
while too small sizes less than 20 pixels are not suggested 
for its bad calculation stability. In fact, the detection per- 
formance, in case studies, is almost always improved with 
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the increasing size of the reference windows; however, the 
improvement is not obvious, but computational cost is sig- 
nificant when the size is over 30 pixels. That is, the 40 pixels 
size is enough for a satisfying detection performance. Any 
window larger than this size is meaningless for performance 
improvement but can contribute a lot to computational cost. 

In addition, for overcoming the limitation of the LDDG 
algorithm itself in distinguishing between bright cells from 
dark cells, we have adopted a simple and effective gray- 
level discrimination rule. The empirical formula for adaptive 
discrimination threshold J/h, which is based on the average 
intensity J, of the inputted line, is assigned for ship cells as 


100, Iq < 50 
50 + Ia, 50 < I, < 100 
In ; (3) 
Ia + (200 — Ig) /2, 100 < Ig < 200 
Ta Ia > 200 


2.3 Analysis of computational complexity 


To analyze the computational complexity, the typical two- 
parameter CFAR method [2], as the widely used ship detec- 
tion algorithm on which several current operational systems 
rely, is selected for comparison. Based on the hypothesis 
that the clutter statistics obey Gaussian distribution, the two- 
parameter CFAR method involves three moving windows: 
target window, guard window and background window. The 
detection criterion can be expressed by 


Xt — Ue | > Tcrar Target 
Oc < Tcrar Cluster 
Target 


(for real implementation), 


(4) 


where X; is the intensity of the test cell, yy is the intensity 
mean in the target window, Tc fr, ris a constant that controls 
the false alarm rate. uc and o, represent the background mean 
and the background standard deviation, respectively, which 
can be computed by the following equations 


Cluster 


Ht — He | > TCFAR 
< TCFAR 


Oc 


1 
Me= 5 2*6) (5) 
© i, jEQe 
JE C 
1 na 5m2 
Oc = a >> (x (i, j) — Xe) ; (6) 
© i jene 


where, x(i, j) is the intensity of cell at (i, j). Qe describes 
the statistical region (background region) of ue and oc, and 
Nc is the number of cells in the region. 

Given that the background window size for clutter statis- 
tics is Ng for an N by N cells image, the minimal number of 
cells for clutter statistics is 4(Ng — 1) for the CFAR method 
and Nz for our S-LGGD method. Then the corresponding 
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Table 1 Computational complexity of the two-parameter CFAR method 
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For each cell 


Total for whole image 


hi Oc Comparison 

Addition 4(Ng — 1)- 1 4(Ng — 1) + 4(Ng — 1)-1 1 N?(12Np — 13) 
Multiplication 1 4(Ng— 1)+2 1 4Ng N? 
Table 2 Computational complexity of the S-LGGD method: LGGD detection step 

LGGD detection for each cell Total for whole image 

Comparison Ai exp(Ai) G Gi 
Addition Ng-—1 0 ~ Ng/2 0 0 0~ Ng/2 (Np — 1)N2 ~ (2Nz — 1)N? 
Multiplication 0 0 0 ~ Ng(Ng +2)/8 1 0 N? ~ (N2 +2Npg + 8)N?7/8 


Table 3 Computational complexity of the S-LGGD method: discrimination and confirmation steps 


Gray-level discrimination for each line 


Total Iavg Comparison 
Addition N 0 N 
Multiplication 0 1 0 


Targets confirmation for each line Total for whole image 


Comparison 
N 3N? 
0 N 


complexities can be computed as listed in Tables 1 and 2, 3. 
In Table 2, all the theoretical maximums are computed under 
the assumption that there are at most half of cells in the back- 
ground window in correlation with the test cell (in practice, 
the relevant cells are far less than the assumed number in 
general). 

Therefore, the numbers of addition and multiplication 
operations are (12Ng — 13)N 2 and 4NgN?, respectively, 
for the CFAR method, while the numbers are (Ng +2)N? ~ 
2(Ng+1)N? and N(N +1) ~ 1/8(N3 +2Ng +8)N? +N, 
respectively, for the S-LGGD method. Thus, the ratios of 
them are as follows: 


N? (12Np — 13) N? (12Np — 13) 


R = 
Sda N? (Ng — 1) N2 QNg — 1) 
~12~6, a 
M 4Ng N? 4Ng N? 
Mul = 
~ NA+ 1 N2(iN2Z+1Ng +1) +N 
32NpB 
x (Ng +4) ~ ——— 2 — ~ (Np +4 
ean? N2 +2Ng +8 EEA 
32 
7 8 
Np 42 (8) 


It is easy to see that the computational complexity of the 
CFAR method is obviously higher than that of the S-LGGD 
method when Ng = 30 is taken into account, even if there 
are always half of cells in background window in correla- 
tion with the test cell for the latter. For any size larger than 
30 pixels (corresponding to 40 pixels window size due to 
Np = 2N = 30 and 2L = 10), further comparison is mean- 
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ingless, because the performance improvement is insignifi- 
cant, but the computational costis significant for the S-LGGD 
method as described in Sect. 2.2. In fact, the computational 
complexity of the S-LGGD method may be always lower 
than that of the CFAR method, since the relevant cells are 
always far less than half of cells in background window, 
which decreases the complexity greatly. 


3 Experimental and results 


To verify the effectiveness of the method described, the S- 
LGGD detector (40 pixels sliding window) has been tested 
over a series of real SAR images, including air- and space- 
borne SAR images. Due to the lack of ground truth data, a 
precise cross-check could not be performed. The validation 
of the detection results was based on the visual inspection. 
For comparison, the two-parameter CFAR method is used, 
and the corresponding results at a false alarm rate of 0.05 % 
are presented. 


3.1 Capability test 


A situation of multiple targets in slightly heterogeneous back- 
ground, as shown in Fig. 2a, is firstly used for the capability 
investigation. The C-band ERS-2 SAR image (400 x 344 
pixels, 12.5m pixel size) shows a busy port containing 49 
ships (bright features) in heterogeneous regions. It should be 
a difficult task since different scales of ships are included 


SIViP (2015) 9:1975-1982 


1979 


Fig. 2 Ship detection for multiple targets situation. a Original SAR 
image; b CFAR result (target window: 3 x 3 pixels, guard window: 
19 x 19 pixels); c CFAR result (target window: 5 x 5 pixels, guard 


in the image, even some small, faint and hidden ships. 
Figure 2b-f presents the detection results by the CFAR 
method and the S-LGGD method respectively. 

Due to different scales of ships in the image (from about 3 
cells to over 20 cells), the choices of detection windows are 
difficult for the CFAR method. For objectivity’s sake, various 
CFAR windows have been tested on the image for obtaining 
an optimal detection result. Figure 2b—e shows the four better 
results. As can be seen, the small target window is helpful for 
detection of small ships, while the large window is helpful 
for detection of large ships (Fig. 2b, c). In addition, with the 
increasing guard window size (from 19 to 31) that is closer 
to the two times the most ship sizes, the detected false alarms 
gradually decrease. 

To sum up, despite some faint ships hiding in the sea clut- 
ter background, the proposed detector highlights all the 49 
ships except for the only one false alarm (as encircled by 
white circle in Fig. 2f). In contrast, the CFAR method causes 
more false alarms even for the optimal result (Fig. 2e). It not 
only causes four false alarms (as encircled by white circle in 
Fig. 2e) but also misses two faint ships numbered as 1 and 2 
in Fig. 2a. Moreover, it is worth noting that a nearly perfect 
shape detection is achieved in comparison with the CFAR 
method, which is vital to the further estimations of the ship 
parameters, including length, width and orientation. 


window: 19 x 19 pixels); d CFAR result (target window: 5 x 5 pixels, 
guard window: 25 x 25 pixels); e CFAR result (target window: 5 x 5 
pixels, guard window: 31 x 31 pixels); f S-LGGD result 


Some closely separated ship situations often result in a 
recurring detection loss to some detectors, such as CFAR- 
like detectors, due to one or more interfering targets occupy 
some of the reference window cells [13]. To test the situa- 
tion, an ERS-2 SAR image of Suez Bay in north of the Red 
Sea (Fig. 3a, 720 x 720 pixels) is adopted. Two white rec- 
tangle boxes indicate two noticeable areas A and B in which 
many ships jammed. The corresponding detection results by 
the CFAR method and the S-LGGD method are shown in 
Fig. 3b—e. As can be seen, the CFAR method with the optimal 
window misses two ships in the jamming area A and one ship 
in the jamming area B, whereas, the S-LGGD method finds 
all ships successfully in both test areas of concern, which 
means that the proposed detector has an adequate level of 
immunity to adjacent-target interference. 

Strong noise background situation is always a challenge 
for most detectors, since ships will be mixed with the ocean 
clutters due to strong backscattering echo of them. To test this 
situation, an ALOS HH polarized L-band SAR image (307 
x 317pixels, 10m pixel size) is adopted in this paper. Fig- 
ure 4 shows the detection results. As can be seen, the hidden 
ship can be detected perfectly by both of the CFAR method 
and the S-LGGD method except for one false alarm caused 
by the latter. However, it should be acceptable for a detec- 
tor in pursuit of speed and automation, since false alarms 
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Fig. 3 Ship detection for closely separated target situation. a Original SAR image; b CFAR result of area A; c S-LGGD result of area A; d CFAR 


result of area B; e S-LGGD result of area B. 


are always inevitable for all the existing detectors in such 
situation. 

The statistical results of Figs. 2 and 3 obtained by the two 
methods are separately listed in Tables 4 and 5. Four statisti- 
cal parameters are used to analyze the detection results, and 
figure-of-merit (FoM) [14] is used to assess the detection 
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and lower false alarm rate. It can be defined as follows: 
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Fig. 4 Ship detection for strong noise situation. a Original SAR image; b CFAR result; c S-LGGD result 


Table 4 Statistical result of SAR imagery 1 


Missed number Detected number False alarm number FoM 

Two-parameter CFAR method 47 4 0.922 

The proposed S-LGGD method 49 1 0.980 

Table 5 Statistical result of SAR imagery 2 

Missed number Detected number False alarm number FoM 

Area A Two-parameter CFAR method 2 22 0 0.916 
The proposed S-LGGD method 0 24 0 1.00 

Area B Two-parameter CFAR method 1 15 0 0.938 
The proposed S-LGGD method 0 16 0 1.00 


where N; counts the accurate detected ship number, and N fa 
counts the false alarm number, Ng; counts the real ship target 
number. 


3.2 Speed test 


The speed test was performed on an SAR strip image acquired 
by a Ku-band airborne SAR system, which has a swath width 
of 8,192 pixels and a resolution of 1.5 mx 1.5 m for both range 
and azimuth direction. The computer on which the detec- 
tion is implemented, configured with an Intel Core 2 Quad 
CPU running at 2.53 GHz, is a general purpose computer. 
The tested data processing rate is near 15.6Mbps for the 
CFAR method, but reaches approximately 38 Mbps for our 
S-LGGD method. Compared with some well-known space- 
borne SAR missions, it is roughly four times faster than the 
designed data downlink rate (10 Mbps) for high-resolution 
radar imaging of NASA’s planetary imaging radar missions 
[15] and near the real-time data processing rate of 45 Mbps 
for the SIR-C/X-SAR mission [16]. 


In fact, a lower data downlink rate, compared with space- 
borne SAR systems, is universal for airborne SAR systems. 
In addition, down-sampling of the downlink image is also the 
universal preprocessing for real-time full-screen monitoring 
in practical engineering application. Taking four to one sam- 
pling ratio into account, the equivalent processing rate for 
the S-LGGD detector can be improved to 150 Mbps, which 
is enough to meet the needs of the real-time ship detection for 
most airborne SAR systems, even for some space-borne SAR 
systems if some special purpose processors are available. 


4 Conclusions 


A linear LGGD algorithm for ship detection from SAR 
imagery is proposed in this paper, and the algorithm uses the 
local gray-level gathering characteristics of SAR imagery to 
characterize its spatial intensity distribution, independent of 
any prior knowledge about ships and background observation 
hence is able to execute automatic ship detection. Based on 
this, a fast and effective detector for ship targets, S-LGGD for 
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short, has been developed. The scheme is easy to implement 
and particularly well suited for a fast stream processing of 
SAR images, which suggests that it has the ability to support 
the real-time ship detection applications. 

The S-LGGD detector has been tested on both air- and 
space-borne SAR images. Due to its sensitiveness to the 
intensity difference, the S-LGGD method has strong differ- 
entiation ability for the target edge. When the detection win- 
dow size is larger than 30 pixels, its detection performance 
is robust for the most situations and has little dependence on 
the sliding window size, since the result is always stable. The 
performance test has shown that the detector can highlight the 
cells with unusual gray-level in sliding window and hence is 
sensitive to ship targets, even to some faint targets. Also, the 
detector is immune to the adjacent-target interference in SAR 
imagery thus is particularly suitable for the multiple targets 
situation. The speed test demonstrates that the detection rate 
of the S-LGGD detector can meet the needs of any real-time 
SAR systems with a data downlink rate under 38 Mbps (or 
150 Mbps for 4—1 sampling ratio processing mode). How- 
ever, the detection performance of the proposed detector is 
decreased with the increasing image noise and heterogeneity, 
thus the detector still needs a further improvement for wide 
applications, especially for some strong noise or heteroge- 
neous situations. 
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